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SOLUTION TO THE VOLTERRA INTEGRAL EQUATIONS OF 
THE FIRST KIND WITH PIECEWISE CONTINUOUS 
KERNELS IN CLASS OF SOBOLEV-SCHWARTZ 
DISTRIBUTIONS 

DENIS N. SIDOROVi'2 



. Abstract. Sufficient conditions for existence and uniqueness of the solution 

' of the Volterra integral equations of the first kind with piecewise continuous 

^ |. kernels are derived in framework of Sobolev-Schwartz distribution theory. The 

' asymptotic approximation of the parametric family of generalized solutions is 

constructed. The method for the solution's regular part refinement is proposed 
using the successive approximations method. 



1. Introduction 

Let us define the triangular region D = {s,t]0 < s < t < T} and introduce 
the functions s = ai(t),i = l,n, which are continuous and have continuous 
derivatives for t G (0,T). We suppose aj(0) = 0, < < ■ ■ ■ < a„_i(t) < t 

for t G (0,T), < «'i(0) < ■ ■ • < a^_i(0) < 1, and functions s = ai(t), i = 0,n, 
>■ . tto(t) = 0, an{t) = t, split the region D into the following disjoint sectors Di = 

a^ _ 

i {s,t : < s < ai{t)}, Di = {s,t : < s < ai{t),i = 2,n}, D = [j Di. Let us 

•n . 1 

*^ ' introduce the continuous functions Ki{t, s) defined for t,s E Di, and differentiable 

^ ■ wrt t, i = 1, n. 

^ ! Let us consider the integral operator 



t „ o.(t) 







with piecewise continuous kernels 

r K,{t,s), t,seDi, 
K{t,s)=l (2) 

[ Kn{t,s), t,S e Dn. 

In this paper we deal with the following Volterra integral equation 

t 

j K{t,s)u{s)ds = f{t), <t <T < oo, (3) 
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where function /(t) has a continuous derivative for t e (0, T), /(O) 7^ 0. Equation 
(3) we call the Volterra integral equation (VIE) with piecewise continuous kernel. 
Our objective is to construct the solution of VIE (3) in the space of Sobolev- 
Schwartz distributions [3]. Obviously, VIE (3) does not have classic solutions 
since /(O) ^ 0. 

The differentiation of VIE (3) leads to integral-functional equation and its 
solution in not unique in the general case [16]. That is why study of VIE (3) can 
not be performed using only the classic methods in the Volterra theory [1,6,15,14]. 
In this paper we continue our results on VIE studies [11], [8,9], [17,19]. We 
consider the equation (3) using the elementary results of the theory of integral 
and difference equations, functional analysis, Sobolev-Schwartz distributions and 
theory of functional equations with perturbed argument of neutral type [7] . 

This paper is organized as follows. Section 2 outlines the construction of the 
singular component of the solution and the integral-functional equation for the 
regular component of the solution is derived. In section 3 we obtain the sufficient 
conditions for existence and uniqueness of solution of VIE (3) in the following form 
u(t) = a6(t) + x(t), were 6(t) is Dirac delta function, x(t) is regular continuous 
function. Such solutions satisfy to the equation (3) in sense of Sobolev-Schwartz 
distributions [3]. To the best of our knowledge, similar studies on VIE (3) have 
not yet been reported in literature. In section 3 we construct the regular part 
of the solution using the the "step method" [13] from the theory of functional 
equations and successive approximations method. In sections 4 and 5 we address 
the most interesting case when VIE (3) has family of solutions depending on 
free parameters. The method for construction of asymptotic approximations of 
parametric solutions is proposed and iterative refinement method is constructed. 
It is to be noted that known method of A.O. Gelfond (readers may refer to [4], 
p. 338) of solution of difference equations is employed. 



2. Definition of the singular component of the solution 

Let us extend f(t) on negative semiaxis with zero and differentiation of VIE 
(3) yields the following equivalent equation 

Fiu)"^ K^{t,t)u{t) + Y,<it)\K^{t,a,{t))- (4) 

n 

-K,+i(t,a,,(t))|u(ai(t)) + ^ j Kl'\t,s)u{s)ds = f^'\t) + f{0)5{t), 

where ao = 0, a„(t) = t. Le us assume i^i(0,0) 7^ 0, Kn(t,t) 7^ 0, for t G [0,T]. 
Let us introduce the following functional operator 

n-l 
i=l 
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and integral operator /Cm = ^ J K~^{t,t)Kl ' {t, s)u{s)ds. 

* = lQi_l(i) 

Then equation (4) can be reduced to the following equation 

u{t) + Au + Ku = K-\t, + K-\0, 0)/(0)5(t). (5) 

Let us search for a solution to VIE (3.2) of the form u{t) = a6{t) + x{t), where a 
is constant, x{t) e C(o,t)- 

It is easy to verify the following identities: 

Ql(t) 



ai{t) 

J ^^^(^M^ = o 



for i = 2,n. Indeed, the first identity holds because > 0, ^^^^^6{s) = 

ai{t) 

^^^^^6{s), J S{s)ds = 9{ai{t)) = 1 for t > 0, were 9 is Heaviside function. The 



second identity is becomes trivial if we notice that for i = 2,n supp 5{s) (iDi = 0, 
J 5{s)ds = 9{ai{t)) - 9{ai^i{t)) = 0, since < ai{t) < a2{t) < ■ ■ ■ < an{t) = 

t. Let us also recall the identity 6{ai{t)) = (refer, e.g. to [3], p. 34). Let us 

take into account the outlined identities and substitution u = a6{t)+x(t) leads the 
equation (5) to the following equation ^"^0, 0)Ki(0, 0)a6{t)+K-^{t, t)^^^^a+ 
x{t) + Ax + Kx = K-\t,t)f^^\t) + K-^{0,0)f{Q)6{t). Equating the last equation 
coefficients of S{t) result a = jJ^^^q^- It is remains to determine the regular part 
from the equation 

x{t) + Ax + Kx = J{t), (6) 

where J{t) = K-\t,t)!^f^\t) - ^^%r^ ;^f|§oy | • It is to be noted that due to the 
operator equality 

Kn{t,t){I + A + K)x = F{x) 
the equation (6) can be written as follows 

Fix) = fit) - (7) 
dt i^i(0,0)- 
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3. Sufficient conditions for the existence of a unique generalized 

solution 

Since 7^1(0,0) ^ then homogeneous equation (4) has only trivial solution of 
singular functions 

m 

def 



Using 





with support at the origin. Therefore, the existence and uniqueness of generalized 
solutions of the equation (4) 

U{t) = Using + X{t), 

x{t) e C(o,T) is equivalent to proving the existence of unique solution of equation 
(6) in C(o,T)- Let us introduce the function 

n-l 



\Ait)\ 



dcf 



i=l 



Let the following condition be fulfilled 



(A) |A(0)| < 1, sup \K-\t, t)K {t, s)\ <c< oo. 

0<s<t<T 

Condition (A) is fulfilled if a^/^\o) are sufficiently small. Here and below the 

n n 

kernel K{t,s) in [jDi is defined as (2). It's derivative wrt t for t, s G IJ-Dj is 
1 1 

defined as follows: 



s 



Theorem 3.1. (Sufficient conditions for existence and uniqueness of generalized 
solutions). 

Let the condition (A) be fulfilled, Ki{t,s) in (2) are continuous functions, are 
they have continuous derivatives wrt t, function f{t) has continuous derivative, 
/(O) 7^ 0. Let i^i(0,0) 7^ 0. Then equation (3) has the unique solution 

where x(t) G C(o,t)- the same time we can find x{t) using the step method 
combined with successive approximations. 

Proof. Since the singular part of the solution is defined then let us consider the 
equation (6) satisfied by the regular component x(t). 

Let us fix g < 1 and select hi > such as sup \A(t)\ = q < 1. Due to 

0<t<hi 

the condition (A) such a variable hi > exists. Let < h < min{/ii,i^}, 
where variable c is defined in condition (A). Let us divide the interval [0,T] into 
subintervals 

[0, h], [h, h + eh], [h + eh,h + 2eh], .... (8) 
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We denote by Xo(t) the restriction of the solution x{t) into [0,h], and by Xm(t) 
we denote it's restriction into subintervals 

Im = [(1 + (m - l)e)h, (1 + me)h], m = 1, 2, . . . . 

Let us select e from (0, 1] such as for t E Im "perturbed" arguments ai(t) G 

m-l 

U 4, z = l,n- 1. If < al \t) < Y^fort G [0,T), i = l,n- 1, then the above 
fc=i ^ 

inclusion holds in the interval [0, T). This inclusion makes it possible to apply the 
well-known in the theory of functional differential equations "step method" . The 
readers may refer to [13], p. 199. 

Let us construct the sequence {xo(t)} : 

x^(^) = -Ax^^-i^xrl+7(^), 

x'o{t)=J{t),te[o,h]. 

to define xo{t) G C[o,fc] 

Due to the selection of h we have an estimate \ \A + K\\c{Ci^o ^^^c^o h)) ^ 
Therefore for t G [0, h] exist unique solution Xo{t) of equation (6). The sequence 
(t) uniformly converge to the solution. We continue the process of constructing 

the desired solution for t > h, i.e. on the intervals In, n = 1,2, ... . For sake of 

clarity let e = 1 in (8). 

Once we get the element xo{t) G C[o,/i] computed we will look for an element 

Xi{t) in the space C{h.2h)- We will find xi{t) from the Volterra integral equation 

of the 2nd kind 

t h 

x{t) + J K;;\t,t)K[{t,s)x{s) ds = J{t)-Axo-j K-\t,t)K[{t,s)xo{s) ds 

h 

using the successive approximations, with XQ{h) = xi{h). 
Let us introduce the continuous function 

^ ~ 1 xi{t), h<t<2h, 

which is the reduction of continuous solution x{t) on to [0, 2h]. Then we can find 
element X2(t) G C(2h,3h) using the successive approximations from the Volterra 
integral equation of the 2nd kind 

t 



x{t) + j K-^{t,t)K[{t,s)x{s)ds 



2h 

2h 



f{t) - Axi- j K-^{t, t)K[{t, s)xi{s) ds. 



Continuing this process for TV steps (A^ > |-) 

We will construct the desired solution x{t) G C(o,t) of VIE (3) using this process 
for N steps. 

This completes the proof of the theorem. □ 
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4. Construction of an asymptotic approximation x(t) of the 

REGULAR PART OF PARAMETRIC FAMILY OF THE DESIRED SOLUTION 

Let us consider the equation (7) which is satisfied by the regular part of gen- 
erahzed solution. Let the following condition be fulfilled 

N 

(B) Exist polynomials Vi = Yl Ki^^f^s^, i = l,n, 

N N 



/ W = E f'^'t"^ Ciiif) = E cuvt^'^i = 1,^-1, where < an < a2i < 

■ ■ ■ < a„_i 1 < 1, such as for t — )■ +0, s +0 we have the following estimates 
\K,it,s) -V.it, s)\ = oat + sf+'), t = \f{t)-rit)\ = 
= C(t^+i),^ = l,n-L 

Expansion in powers of t, s which are presented in the condition (B) we will 
call as "Taylor polynomials" of the corresponding functions. Let us introduce the 
function 

BU) = J^„(0,0) + ^(a:(0))i+^(i^,(0,0) - J^,+i(0,0)), 

which depends on integer argument j, j G NUO. Function B{j) which corresponds 
to the main "functional" part of the equation (7) is called as characteristic func- 
tion of equation (7). Let us consider the construction of an asymptotic solution 
of equation (7). 

In contrast to section 3, in sections 4 and 5 it is not supposed that homoge- 
neous equation for equation (3) has only trivial solution. Therefore the solution 
of integral- functional equation (7) can be non unique. Let us follow paper [5] and 
search for the asymptotic approximation of a particular solution of the inhomo- 
geneous equation (7) as following polynomial 

N 

x{t) = Y,Xj{\nt)t^. (10) 

j=0 

Let us demonstrate that coefficients Xj depends on Int and free parameters in 
general irregular case. This is consistent with the possibility of the existence of 
nontrivial solutions of the homogeneous equation. 

For computation of the coefficients Xj we consider regular and irregular cases. 

Definition 4.1. Point j* is called regular point of characteristic function B{j), 
if B{j*) 7^ and irregular point otherwise. 

4.1. The regular case: characteristic function B{j) ^ for j G (0, 1, . . . , A^), 
where is sufficiently large. In this case, the coefficients of Xj will be con- 
stant, i.e. independent on Int. Indeed, lets substitute the expansion (10) into 
the equation (7). Using the method of undetermined coefficients and taking into 
account the conditions (B), lead to the recursive sequence of the systems of linear 
algebraic equations wrt Xj : 

/(O) /(O) dK,{t,Q) 



B{0)xo = f (0) 



iri(0,0) i^i(0,0) dt 



(11) 
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=M,(xo,...,x,_i), J = 1,...,N. (12) 

Mj are expressed in terms of solutions Xq, . . . ,Xj_i of previous equations and 
coefficients of the Taylor polynomials from the condition (B). 

Since in the regular case B{j) ^ the coefficients XQ,...,Xfq can be uniquely 
determined and the asymptotics (10) can be constructed by this means. 

4.2. Irregular case: characteristic function -B(j) in (0, 1, . . . , A^) has zeros. 

Let us demonstrate that in irregular case the coefficients Xj are polynomials in 
powers of Int and depends upon arbitrary constants. The order of polynomials 
and the number of arbitrary constants are related to the multiplicities of integer 
solutions of the equation -B(j) = 0. 

Indeed, since the coefficient xq in the irregular case can depend on Int, then 
based on the method of undetermined coefficients x^ can be found as the solution 
of the difference equation 

n-l 

ir„(0, ^)x^{z) + «^(0)(^^(0, 0) - i^,+i(0, 0))xo(^ + a,) = (13) 
^^,(Q) /(O) K^{t,^) 



Ki(0,0) dt 

where Oj = lnQ;'(0), z = \nt. There are three possible cases here: 
1st case. (5(0) ^ 0). 

In this case the coefficient xq does not depend on z and can be determined 
uniquely from the equation (11). 
2nd case. (5(0) = 0). 

Let j = is simple zero of the function B{j), i.e. -B(O) = 0, -B'(O) ^ 0. Then the 
coefficient Xq{z) we can find from the difference equation (13) as linear function 

Xo{z) = XoiZ + Xo2. (14) 

Lets substitute (14) into (13). Thus for determination of the coefficients xoi,xo2 
we obtain two equations as follows: 

5(0)xoi = 0, (15) 

(16) 



t=0 



where -B(O) = 0, -B''^^(0) ^ 0. Hence the coefficient xo{z) is linear wrt z and 
depends on the arbitrary constant. So, it the 2nd case 



^»<^' = (^"(«' - tm^r^Jwwr ^ 

where c is arbitrary constant. 

3rd case. Let j = be root of the equation B{j) = with order of multiplicity 
of A; + l, i.e. 5(0) = 5'(0) = ...5(*^)(0) = 0, 5('=+i)(0) ^ 0, A; > 1. Solution xq (2) 
of the difference equation (12) we search in the form of a polynomial 

Xo{z) = Xoiz''^^ + X02Z'' H h Xok+lZ + Xok+2- (17) 



D.N.SIDOROV 



Let us substitute the polynomial (17) into the equation (13) and take into account 
the identity 



n-l 



B{j) = J](a:(0))i+X'(^.(0,0) -ir,+i(0,0)), 



where = lnQ;^(0). Next lets equate the coefficients of powers 



z ^ z ^ . . . ^ z 



to zero. Finally we get recurrent sequence of linear algebraic equations wrt 

{ 5(0)xoi=0, 

i?(0)xo2 + i?(^)(0)f Voi = 0, 



k 

B (0)xo.+i + ^^'^0) ( A: + 1 - / ) + ^^'"'^ ( A: + 1 - / ' 

■■■ + i?(^no)f ^t!7^ + M^o. = o,/ = i,...,fc 



S(0)xofc+2 + 5('=+i)(0)xoi + 5«(0)xo2 + . . . 5(i)(0)xofc+i 

/(O) 9A'i(t,0) 



(19) 



In our case 5(0) = 5'(0) 
mial (17) we let xqi 



B('=+i){0) 



iri(0,0) 9t 
= 5(^')(0) = 0, 5(^+1) (0) 7^ 0. Hence in polyno- 
f'^^) ~ x((o 0) ) • Equations of the system 



(18) become identities i?(0)xoj = 0, j = 1, /c + 1, -B(O) = 0. Hence the coefficients 
X02-I ■ ■ ■ ,xok+2 of the polynomial (17) remains arbitrary constants. 

Next, let's employ the method of undetermined coefficients and take into ac- 
count the identity 



f t^ ln'^trft = t^+^^(-l 



(j + 1) 



s+l 



By this means we construct the difference equations for determination of the 
coefficient Xi{z) {z = Int) and next coefficients of the asymptotics (10). Indeed, 



L(x] 



def 



x=xo(z)+xi (z)t 



n-l 



i^„(0, O)xi(^) + 5^(a:(0))2(J^,(0, 0)- 



(20) 



i=l 



-Ki+i{0, 0))xi(z + a,) + Pi{xo{z)) 



t + r{t), r{t)=o{t). 



Here Pi{xq{z)) is the polynomial of z, which degree is equal to the multiplicity of 
solution j = of the equation i?(j) = as have been proved. From the relation 
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(20) due to r{t) = o{t) for t — >■ it follows that coefficient Xi{z) have to satisfy 
the difference equation 

n-l 

K„(0, 0)x,{z) + ^(a'(0))2(i^,(0, 0) - K,+,{0, 0))xi{z + a,)+ (21) 



i=l 



+Pi(xo(z)) = 0. 

If -B(l) 7^ 0, then the equation (21) has the solution Xi{z) as the same degree 
polynomial as multiplicity order of the solution j = of the equation B{j) = 0. 
If j = 1 is also the solution of the equation B{j) = the solution xi{z) can 
be constructed as polynomial of the power + ki, where ko and ki are the 
corresponding multiplicities of the solutions j = and j = 1 of equation B{j) = 0. 
Coefficient Xi{z) depends on k^ + ki arbitrary constants. 
Let us introduce the condition 

(C) Let the equation B{j) = in the array (0,1,..., A^) has the solutions 
ji, . . . ,ju of multiphcities ki, i = l,u. 

Then, in a similar way we can calculate the remaining coefficients ^2(2;), . . . , Xis[{z) 
of the asymptotic approximation x{t) of solution of equation (7) from the follow- 
ing sequence of difference equation 

ir„(0, Q)x,{z) + ^(a'(0))i+^^ (i^,(0, 0) - i^,+i(0, Q))x,{z + a,) + 

i=l 



+V,{xq{z), . . . , Xj.i{z))) = 0, J = 2, N. 
Thus we have the following lemma 



Lemma 4.2. Let conditions (B) and (C) he fulfilled. Then exists the function 

N 

x{t) = '^Xi{lnt)f, such as for t — )■ +0 the residual solution of equation (7) 

1=0 

satisfies the estimate 



F(x(t))-/(i)(t)+ir«(t,0)-/^°^ 



A'i(0,0) 



o{t 



The coefficients Xj(lnt) are polynomials of Int. The degrees of these polynomials 
are increasing and do not exceed the sum of the multiplicities of ^ kj of solutions 

j 

of the equation B{j) = from the array (0,1, dots,i). The coefficients Xi{lnt) 
depend on Ylj=o arbitrary constants. 

i 

Remark 4.3. If B{j) 7^ 0, then in the sum ^ kj we zero the corresponding kj. 

j=0 
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5. An existence theorem for continuous parametric solutions 

families 



Since < a'-{0) < 1, ^^(O) = 0, z = l,n — 1, then for any < e < 1 exists 
T' G (0, T] such as the following estimates are fulfilled 

ma x < 

i=l,n~l,te[0,T'] 

sup < e. 

i=l,n-l,te{0,T'] ^ 

Let us introduce the condition 

(D) Let Kn(t,t) 7^ for t e [0,T'] and A^* is chosen so large that the following 
equality 

sup e^*\A{t)\ <q<l 

te{0,T') 

is fulfilled, where function A{t) is defined the Section 3 with formula (*). 

Lemma 5.1. Let the condition (D) he fulfilled. Let in C(o,t') class of continuous 
functions fort e (0,T'] which have the limit (which could he infinite) fort — +0 
exists an element x{t) such as for t — )■ +0 error of the solution of equation (7) 
satisfy the estimate 

\F{x{t)) - fit) + ^((^'0)-^^| = o(t^), 

where N > N*. Then equation (7) in C(o,t') has the solution 

x{t) = x{t) +t^v{t), (22) 

where v{t) is uniquely determined hy successive approximations. 

Proof. Substitution of (22) in equation (7) gives us the following integral-functional 
equation for determination of the function v{t) 

v{t) + K^it,t)I^J2<^^^{^y {K,{t,a^{t))- (23) 
-K,^i{t,a,{t))]viaiit)) + J2 J Kl'\t,^)iA] <s)ds 

Let us introduce the linear operators 

Mu K-\t, t) 1^ a[{t) I A,(t, am- 

-Ki+i{t,ai{t)) \v{ai{t)), 



VOLTERRA EQUATIONS WITH PIECEWISE CONTINUOUS KERNELS 11 
n '^f 

Kv'^Y. / K-\t,t)K['\t,s){s/tf\is)ds. 

Then equation (23) can be presented as following operator equation 

u+{M + K)u = 7(t), 

where ^{t) is the right hand side of the equation (23). This function is continuous 
due to the condition of the Lemma 2. Let us introduce the Banach space X of 
continuous functions v(t) with norm 

\\v\\i = max e~'-^\v(t)\, I > 0. 

0<t<T' 

Then due to the inequalities sup < £^ < 1 and due to the condition (D) 

t€{0,T'] 

for V/ > norm of a linear function of the operator M satisfies 

\\M\\cix^x) <q<l. 

In addition, for the integral operator K for sufficiently large / the following esti- 
mate is correct 

\\K\\c{x^x) < gi < 1 - g- 

For sufficiently large / > this implies that 

\\M + KWax^x) < 1, 

i.e. the linear operator M + is a contraction operator in the space X. Hence 
the sequence {vn} converge where Vn = — (M + K)vn-i + 7(^), = 7(t). This 
completes the proof of the theorem. □ 

Theorem 5.2. Let the following conditions be fulfilled (B), (C), (D),f {0)^0, 
Ki{0, 0) 7^ 0. Then equation (3) for < t < T' < T has the solution 

x{t) = j^^^6{t)+x{t)+t'''v{t), 

V 

which depends on ^ ki arbitrary constants, where ki are determined in the con- 

i=l 

dition (C). Function x is constructed in the form of (10), then v(t) is uniquely 
determined with successive approximations. And we have the following asymptotic 



estimate 



0{t^ ) fort^ +0. 



Proof. Based on the Lemma 1 because of the conditions of the theorem is possible 
to construct an asymptotic approximation of the regular part x(t) of the solution 
in the form of the following log-power polynomial: 

N 

^Xi(lnt)f. 

In this case, by construction, the coefficients Xj(lnt) depend on the certain number 
of arbitrary constants. Due to the lemma 2 the substitution x{t) = x{t) + t^''u{t) 
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enable the construction of the continuous function u{t) using the successive ap- 
proximations method. □ 

The solution constructed on [0, T'] can be extended on the whole interval [0, T], 
based on known "step method" [13, c.l99]. 

In simple cases one can use the solution of the equivalent equation (4) in order 
to construct the solution of the integral equation (3) in closed form. 

Example 5.3. 

t/2 t 

x{s)ds + 2 j x{s)ds = 2 + t, t > 0. 

t/2 

lYamTilfi Viae fVifi fnlln^srincr fnrm — , , 

2'^y2' 



An equivalent equation (4) in this example has the following form —^x{^) + 



2x{t) = 2S{t) + 1. The desired solution is as follows x{t) = 26{t) + 2/3. 
Example 5.4. 

t/2 t 



/.(.M.-/.(.M.^i + M>o. 



t/2 

Here an equivalent equation is as follows = 6(t) + l. It has c-parametric 

^ ~ In2' 



family of generahzed solutions x{t) = S{t) + c — jsl, c is constant. 
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